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THEORETICAL STUDY OF THE TUNNEL-BOUNDARY LIFT INTERFERENCE DUE TO SLOTTED 
WALLS IN THE PRESENCE OF THE TRAILING- VORTEX SYSTEM OF A LIFTING MODEL ^ 

By Clarence W. Matthews 


SUMMARY 

TJte equations ■presented in this report give the irderjerence on 
the irailing-vortex system of a uniformly loaded finite-span 
vnng in a circular tunnel containing partly open and partly 
closed walls, with special reference to symmetrical arrangements 
of the open and closed portions. Methods are given for extend- 
ing the equations to include tunnel shapes other than circular. 
The rectangular tunnel is used to demonstrate these methods. 
The equations are also extended to nonuniformly loaded wings. 

An analysis of the equations for certain configurations has 
shown that: (1) only a small percentage of slot opening is 
required to give zero interference conditions if the tunnel con- 
tains four or more slots; (S) in the configurations studied, the 
ratio between the slotted-tunnd interference and the dosedr- 
tunnel inteiference at the center of the tunnel is approximately 
constant for various model spans; and (S) tunnds containing 
an odd number of slots or nonsymmetrical slot arrangements 
cause an additional rolling moment or a crossflow on the wings, 
or both. 

INTRODUCTION 

In a study of solid-blockage interference (see ref. 1), it has 
been sbovm that tunnels containing mixed boundaries, that 
is, partly open and partly closed walls, wiU. eliminate or 
greatly reduce such interference. Since the slotted tmmel 
configuration required to eliminate solid blockage can not 
eliminate lift interference, it is necessary to study the inter- 
ference on the trailing vortices of a lifting model in order to 
make the necessary corrections to the lift characteristics of a 
model. 

The problem of one or two slots has been treated by various 
authors (see refs. 2 to 4). The case of more than two slots 
has also been treated in references 6 and 6. Reference 5 
treats only small wings in circular tunnels, and reference 6 
treats only the case for a large number of evenly spaced slots. 

The purpose of this investigation is to present equations 
which express the tunnel-wall interference due to mixed open 
and closed boundaries in the presence of the trailing-vortex 
system of a finite-span lifting model at subsonic velocities. 
Special attention is given to test sections in which the slots are 
symmetrically located with respect to both axes. 

Various extensions of the theory have been made and 
follow in a general fashion the methods of reference 5. These 
extensions include the effects of wing span, slot configurations. 


interference at points near the center, nonuniform loading, 
and methods for calculating the interference in tunnels of 
other than circular cross section. 

Numerical calculations of the interference characteristics 
of several symmetrical cases are presented and are used to 
show the properties of the interference of circular t unnel a 
containing 1, 2, 4, 8, and 12 slots symmetrically located with 
respect to the x- and y-axes and a square tunnel containing 
8 slots symmetrically located in the top and bottom walls. 

SYMBOLS 

An, Bn,\ 

^n, / constants 
0^1 fin / 

b span of wing 

s half-span of wing (see fig. 1) 

C tmmel cross-sectional area 

Cl lift coefficient of model 


s 



AOd drag increment due to interference 
D tunnel diameter (see fig. 4) 

ds tangential line element 

k quality factor, the ratio at any given point of the 

interference of a slotted tunnel to the inter- 
ference of a closed tunnel with the same cross 
section 

I half-span length of wing in f-plane (fig. 1) 

m number of slots or panels in tunnel 

9, 3, 2 integer denoting indices of summations and mul- 

tiple products 

p special number defined in expression (2) 

r radial distance of point from coordinate center 

JR radius of circular tunnel 

S model or wing area 

u a;-component of velocity 

V y-component of velocity 

V magnitude of velocity 

V velocity vector 

X, y coordinates of rectangular system (fig. 1) 
z complex variable, x-\-iy 

^ complex velocity, u—iv 


1 Supersedes recently declassided NAOA BM L63A26, 1953. 
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r circulation about a point, positive in counter- 

clockwise direction 

7 function of 0, and b (see eq. (32)) 

8 proportionality factor used in equation for ex 

ej interference correction angle to be added to meas- 

SSO 

ured angle of attack, (see ref. 7) 

<r strength, of a source 

f complex variable in transformed plane 

6 angular coordinate of polar system (fig. 1) 

dg polar angle of slot edge 

complex function of z (z=e*’) 

X parameter defined in equation (106) (also see para- 

graph 14.8, ref. 8) 
perturbation potential 

by derivative of potential in direction normal to a 

b/i given line 

Subscripts; 

0 closed tunnel 

/ interference condition 

N nommiformly loaded 

U uniformly loaded 

THEORY OF LIFT-VORTEX INTERFERENCaE IN SLOTTED 

TUNNELS 

GENERAI. ANAITSIS 

Theoretical boundary conditions of flow about trailing 
vortices in a tunnel containing mixed open and closed 
boundaries. — ^The equations for the interference on the lift 
of a model due to mixed open and closed tunnel walls in the 
presence of a trafling-vortex system may be obtained by 
considering the same two-dimensional approximation of the 
flow fleld that is used in reference 7. The conditions of this 
two-dimensional approximation may be briefly stated as: 
(1) the tunnel and its boundaries extend from a point an 
infinite distance upstream of the model to a point an infinite 
distance downstream of the model; (2) the velocities induced 
by the trailing vortices in the cross-sectional plane located at 
the model are one-half of those induced in a far-downstream 
cross-sectional plane; (3) the induced-velocity flow field 
in the far-downstream section may be treated with two- 
dimensional methods; (4) the boimdary condition which must 
be satisfied at a solid portion of a timnel wall is that the flow 
must be tangential to it, or 



a) 


(5) the condition which must be satisfled at an open portion 
of a tunnel wall is that the potential over that portion must 
be constant or the flow must be normal to it; (6) no singulari- 
ties other than the trailing vortices can exbt within the 
boundaries of the timnel; and (7) the constant potential in 
every slot must be equal to that in every other slot. The 
final condition is required because the pressure is the same at 
every slot and is shown by the following considerations. 
Since the pressure is constant over the entire region outside 
the tunnel there will be no pressure differences between the 
slots at a point far upstream, and hence there can be no flow 


between the slots due to external influences. Also, since this 
point is too far from the model to be influenced by it, there 
can be no flow due to the model. As there is no flow between 
the slots, no potential gradient can exist between the slots; 
therefore, the potentials in all the slots must be the same. 
It therefore follows by condition (5) that, in the far-down- 
stream position, the constant potential in every slot is equal 
to that in every other slot. 

Coordinate system used. — ^The coordinate system used in 
reference 8 is also used throughout this report. In this 
system, it is expedient, in order to avoid confusion with the 
customary notations, x+iy, to use as axes of reference the 
following system: x-axis to the right, z-axis downstream, 
and j/-axis to form a right-hand system. Since the z-axis 
is not used in the calculations, it is possible to use the comple.x 
coordinate z=x+iy without confusion. The velocity com- 
ponents are denoted by m in the x-direction and v in the y- 
direction. The symbol w is reserved for the complex 
potential. 

Velocity fields in circular slotted tunnels. — The previously 
stated boundary conditions may be satisfied by using 
complex velocity functions rather than complex potential 
functions. This is done by selecting a complex velocity 
function that has singularities at the wing tips and a flow 
direction either normal or tangential to the timnel walls. If 
this function is multiphed by another complex function 
whose value on the tunnel wall changes from all real to all 
imaginary (or the opposite) at each slot edge, then the flow 
will be rotated 90° at those points so that the final flow of the 
product of the two functions wfll be normal to the wall on 
selected portions and will be tangential to the wall on the 
remaining portions. It can, however, be expected that the 
second function will introduce within the tunnel singiflarities 
which are not permitted according to the stated boundary 
conditions, so that a third function which contains all the 
forbidden singularities must be used in such a fashion that it 
will cancel the forbidden singularities of the second function. 

In order to make up the first function, suppose that the 
two singularities at the wing tips with their reciprocal 
singularities are written as l/(z*— s‘)(l — z*s^, whore z is 
the complex variable aj-f-iy and s is the semispan of the lifting 
wing. If this function is examined by letting z=e'®, where 
0 is a polar angle (fig. 1), it wfll be found that the flow may 
be made normal to the walls if the factor z is included in the 
numerator. This flow may also be made tangential to the 
wall by multiplying by the factor i. Thus, the first function 
may be written 


pz 

(z“— 8*)(1— z*s*) 


( 2 ) 


The symbol p can be chosen to be either 1 or i, depending on 
whether normal or tangential flow is required at the walls. 

The second function can be developed by considering the 
square root of a function which is real on the wall and changes 
sign at each slot edge so that the square root of the funetion 
changes from real to imaginary at each slot edge. Such a 




V 


and is in general complex. Examination of this function 
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z-pion'e 

(a) 




(a) Circular tunnel. 

(b) Square t unn el. 

Fioubb 1. — The coordinate system used for the investigation of the lift 
interference due to slotted tunnels. 


shows that it becomes 


V 


cot —cot ^ on the tunnel wall, 


where 0, is the polar angle of a slot edge (see fig. 1). The 
term under the radical changes from positive to negative at 
so that the function changes from real to imaginary at 
the point Q, on the tunnel wall. If a number of these func- 
tions having different Qq's marking the transition from open 
to closed sections of the tuimel "waU are multiplified together, 
it can be seen that the product will change sign at each value 
of Qg so that the function ■will be real on alternate sections. 


This use of cot ^ also suggests the use of other complex trig 


nometric functions such as cos t?, sin t?, and tan ~ The sine . 

and cosine terms 'wiU, on examination, be foimd to introduce 
two slot edges for each value of Qg, rather than one. These 
two slot edges will be found to be located symmetrically with 
respect to the a;-axis for the cosme term and to the y-axis 
for the siue term. Thus, the use of these functions is sug- 
gested for symmetrical slot configurations. 

Since any of these functions may introduce singularities 
within the tunnel, it is necessary to examine them for such 
singularities. The four functions may be written as multiple 


products of the form 
nVcOSt?— = 



Vcost? — — cos 02 • • . Vcost? — COSOq 
with equivalent 'expresssions for 

(3) 

g 

n VsiDt?— sm0, 
0=1 

9 1 s T. 

(4) 

,5.V 

and 

cot --cot ^ 

(^ 

n J 

0-1 V 

4 . 6g 

tan tan 

Jj 2i 

(6) 

When z=e*^ is substituted for i?, the expressions (3), (4), (5), 
and (6) become, respectively. 

n V^+ 1 - 22 COs 0 <, 

00 


g 

n ■y]^—l—2izBui6g 

(8) 

(2_i)-o/s n 
0-1 ' 

J i(z+l)—(z—l) cot^ 

(9) 

(z-f n 

0-1 

^ (z— 1)— i(z-f 1) tan^ 

m 


Examinations of the functions (7) to (10) shows that they 
contain forbidden singularities at the following values of z\ 
the cosine and sine terms at z=0, the cotangent term at 
2 = 1 , and the tangent term at s=— 1. Several functions 
which contain singularities identical to those appearing in 
functions (7) to (10) and which are real on the tunnel are 

yi (aaCOSTOT?-}-6„smTOT?) (11) 

n=0 

which has a singularity at z=0, 

S (a,+il30 cot- (0 (12) 

which has a singularity at z=l, and 

^ (a^+ip^ tan" (13) 

which has a singularity at 2 = — 1. In the expressions (11), 
(12), and (13), a„, bn, and Pn are real constants. It may 

be seen now that if the trigonometric series is di'vided by the 
multiple product which contains the same singularities, the 
forbidden singularities will be canceled out of the final equa- 
tion and an equation ■will be left which has only the desired 
singularities wi thin the tunnel. 

The reason for using the m'ultiple product in the denom- 
inator may be seen by examining the flow about the slot 
edges. Since the flow must turn a sharp comer as it goes 
around the edge of each slot, the velocity must be infinite 
at that point. Since the multiple-product function becomes 
equal to zero at each slot edge, it must be placed in the 
denominator of the final complex velocity function to insme 
the required infinite velocities. 
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Before the final complex velocity function is -written, it 
should be observed that function (2) has a zero at the point 
2=0 -which must also be removed either by including an 
extra term in summation (11) or by multiplying summation 
(12) or (13) by cos sin t?. With this relationship in 
mind, the final complex velocity functions may he -written: 


dw 

dz 


■2+^ 

pz~^ (a„ cos sin n&) 


(2“ — s*) (1—2*8“) n Vcos t?— cos 6g 

ff-i 


a4) 


dw 

dz 



pz (®» cos 7lt?+5a sin Tit?) 

n— 0 


(2“— s“) (l---2“8*) n Vsin t?— sin 0, 

ff-i 


(16) 


summations must be evaluated. This evaluation can be 
effected by using the remaining boundary conditions, which 
require that the potentials m each slot must be equal to 
each other and that only a vortex flow can exist about the 
singularities wi thin the tunnel. 

The potential condition can be evaluated by the foUo-wing 
line integral: 


fCij.l'j) 


( 18 ) 


in which the limits of integration terminate in the slot. One 
path which can be used is ^e streamline which flows along 
the panel. When this streamline is used as the path of 
integration, equation (18) becomes 


«n+l 


V d0=¥>,+i— ¥>,=0 


( 19 ) 


dw 

dz 


pz (oi cos i?+6i sin t?) ^ (o!n+ii3j cot“^0 

(2“— s“) ( 1 — s“8“) n .«/cot^— cot^ 
ff-i V ^ ^ 


ae) 


and an equation similar to equation (16), where tan — is used 

* ^ 
rather than cot-^- 


In the application of these equations, each value of fl, is 
noted to introduce two slot edges into each of equations (14) 
and (15) and only a single slot edge into equation (16). 
Since each panel has two edges, 2 must be equal to the n umb er 
of panels m in equations (14) and (15) and to twice the 
number of panels, or 2m, in equation (16). Also, if equation 
(14) or (15) is applied to tunnels containing an odd number 
of slots, it -will be found that the singularity arising from the 
multiple product -wfll contain a term of the order }i, which 
cannot be removed by the summation. In order to remove 
this singularity, the summation of equation (14) or (15) 
must be re-written as 


ro+l 

s (a„ cos t?+ sin (17) 

This series -wfll, upon examination, be found to contain a 
singiflarity of the order Kj which -will cancel the one-half- 
power term due to the extension of the multiple product 
over an odd number of slots. 

Equations (14) and (15) introduce slot edges at locations 
symmetrical to the x- and y-axes, respectively; thus, it is to 
be expected that the -use of these equations for symmetrical 
cases -wfll simplify the problem. These equations wfll not 
work, however, for nonsymmetrical cases as they -wfll intro- 
duce slots at points where they are not desired; therefore, it 
is necessary for the general case to resort to the more com- 
plicated but completely general equation (16). Equation 
(16) can be used for symmetrical cases; however, the expres- 
sions resulting from its use can be reduced to the forms 
obtained from equations (14) and (15) . 

In order to make the final application of equation (14), 
(15), or (16) to a -wind tunnel, the arbitrary constants of the 


as the potential is the same over all slots. Smce the velocity 
V -wfll in general be complex, two equations exist for each 
panel. 

The strengths of the circulation about a pole and of the 
somce at a pole are fixed by the relation, 

r+i<^=^^d2 (20) 

where r is the circulation due to lift and a is the source 
strength and is equal to zero. Since equation (20) must be 
evaluated about each pole, it gives four equations which 
may be -used to evaluate the constants. 

The set of boundary equations (19) and (20) may be con- 
sidered as a set of simultaneous equations in the -unknown 
constants found in equation (14), (16), or (16) and is used 
to determine the values of these constants. If these con- 
stants are to be uniquely determined, there must be as many 
equations as there are constants. An examination of equa- 
tions (19) and (20) shows that equation (19) gives 2m 
equations and that equation (20) gives four equations, 
making a total of 2m-l-4 equations which may be used to 
determine the imkno-wn constants. The number of constants 
which must be satisfied is determined from an examination 
of equation (16). The paragraph foUo-wing equation (16) 
sho-ws that 2 must equal 2m, so that the number of constants 
must equal 2m-l-4; hence the munber of equations and the 
number of constants are equal, so that all the constants ai'e 
determined and a unique solution is obtained for the problem. 

DEBJVATION OP EQUATIONS 

Slots symmetrically located -with respect to the x- and 
y-axes. — ^In the case of symmetry about the x- and y-axes, 
considerable simplification results. Equation (14) can be 
used in place of equation (16). In equation (14), each value 
of 6, in the multiple products can be used to produce four 
changes of sign, pro-vided these products are -written as 


ml2 

n iJcOS^ 1^—COQ^ Off 
^=1 

(21) 

mt2 

n Vsin“d— sin“0, 
<7=1 

(22) 
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These changes of sign are observed at d,, —dg, Tr—dg, and 
ir-\-6g, or one in each quadrant. Hence, if the values of 6g 
are known in one of the quadrants, the multiple product need 
be extended only over that quadrant, as all the other slots 
win be automatically introduced. If the slots are evenly 
spaced and of equal width, function (21) or (22) can be 
replaced by 


V 


COS'* d— cos * 01 


or 


V . jm 
sm*-2 


m 


d— sm*-^0i 


(23a) 


(23b) 


Expressions (23a) and (23b) can be shown to produce slots 
which are evenly spaced and of equalwidth by the substitution 
of 


where 


m ^ 


^ m 


for d in expressions (23 a) and (23b). 

When this modification is taken into consideration, equa- 
tion (14) then may be written with the substitution of 


s=e*^ as 


dw 


f+i 


2 / 


n-0 \ 




2z» 




2*»— 1 
23* 


) 


1 fn/2 

(3*-s*)( 1-3*«^^57^ n cos 20,+l 


(24) 


which may be reduced to 


dw 


W , , P 

PS ' 

^0 L 


n+l 






--1)] 


. &S) 

( 2 *— 8*)(1— 2 *s*) n -4^—2^ cos 20,-|-l 

17-1 

where the factor p is equal to 1 if slots intersect the aj-axis 
and is equal to i if panels intersect the x-axis. If the slots 
are evenly spaced, and of equal width, the function 


22 " cos TO01-1-1 

can be substituted for the multiple product in equation (24) 
or (25). 

The constants of equation (25) must be evaluated by using 
the boundary condition for the equality of the potential in 
each slot (eq. (19)) and for the circulation and nonexistence 
of sources (eq. (20)). Equation (20) is evaluated by applying 
Cauchy’s integral theorem to equation (25) to obtain the 
required line integral of dwldz. This evaluation gives, about 
the positive pole s. 


r-|-i<r=- 


2irip2^' S [o:^^"'^^(s*“+l)-iM^”'^\s^-l)] 


ftcaO 


2s(l— n Vs^~2s“ cos 20,-+-l 

( 7-1 


(26) 


and about the negative pole —s 


77 >_ 2+1 f H-n+l ^+1 ■) 

2vip22 S I «,(-«) 2 ^ [(-«)^+l]-iU-8)^ [(-s)*’'-l]| 


— r — i(T= 


777/3 ^ 

— 2«(1— 8^ n 2s^ cos 20,-1- 1 


(27) 


The circulation (eq. (20)) is determined by setting F equal to 
the real parts of equations (26) and (27) while the nonexist- 
ence of sources is assured by settmg the imaginary parts of 
equations (26) and (27) equal to zero. The condition on the 
potential (eq. (19)) is satisfied by letting s=e" in equation 
(14), so as to obtain the velocity along the panel streamline, 
which is a convenient streamline to use for the evaluation of 
equation (19). Equation (19) then becomes 


»=x 


f+1 

2D (®i 7 cos nB—ihn sin n6) dd 

wO 

panel cqS 20-ks^ ]I VcOS* 0— COS* dg 

?-l 


(28) 


m/2 


dw 


psr(l— s^) n Vs*— 2s* cos 20J-1-1 

»=i 


dz 


m/2 


where the notation indicates integration over a panel. 
Separate equations are then obtained, one for each panel. 
The constants a, and b„ can now be evaluated by considering 
equations (26), (27), and (28) as a set of simultaneous equa- 
tions in a„ and b„. 

In solving equations (26) to (28) for a„, it is observed that 
a new constant a„ appears; thus, 

a,=a» n Vs*— 2s* cos 2dg+l 


7t2= 


-1 ( 7=1 


(29) 


A corresponding- relation between and a new constant 
also appears. With the substitution of these new con- 
stants an and /3„, equation (25) may be written 


v( 2 *— s*)(l— 2 *s*) n Vz*— 22 * cos 20j-f-l 


sL 


S-n-fl — 

0,2* (2*»+l)— i/3,2* 


n-H 






(30) 
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and equations (26), (27), and (28) become, respectively, 


I t 


-rsL 


n-0 









n+l- 




(31a) 

(31b) 


-X 


y^. (a„ cos nd~i^n sin dd 

n—0 

mf2 


*“**' (1— 2s* cos 20+8*) n V cos* cos* 0, 

g-l 


(31c) 


Examination of equations (31) shows that a number of 
the o,’s and jSn’s are equal to zero. This is shown by first 
considering the nature of the simultaneous equations (31). 
The equating of the real and imaginary parts divides the 
entire set into two separate sets of simultaneous equations, 
one of which involves a^’s only and the other only. "When 
p is chosen, one of the sets will he homogeneous whereas the 
other set will contain two equations with constant terms 
equal to —1. Thus, one set of constants becomes equal to 
zero whereas the other set has values which may he different 
from zero. This remaining set of equations can he reduced 
hy considering the values of the int^rals across any one of 
the panels and its counterpart which is symmetrically located 
with respect to the y-axis. It can he shown that these in- 
tegrals will have the same absolute values and the same sign 
if n is even but opposite signs if n. is odd. The same feature 
is also observed in comparing equations (31a) and (31b); 

that is, if ^—71+1 is odd, then the terms which contain b to 

that power will have opposite signs whereas the remaining 
terms will have the same signs. By adding or subtracting 
equations (31a) and (31b) and the equations for the integrals 
across each of the symmetrically located panels, the set of 
equations in a„ (or in may be reduced to two new sets, 
one of which contains the odd values of n and the other the 
even values of n. Again one of these sets is homogeneous, 
and one equation of the other set has a constant term equal 
to —1 so that all the constants associated with either the 
odd values of n or the even values of n will be equal to zero 
and the other set will have values which can be different 
from zero. 

This analysis shows that one of four sets of constants oay, 
« 27 +ij Pn, mid Pii+i can occur, so that four different symmetri- 
cal slot configurations are suggested. The possible sym- 
metrical configurations which can occur are (case I) panels 
intersecting both axes, (case II) slots intersecting both axes, 
(case m) a panel in tersecting the r-axis and a slot intersect- 
ing the y-axis, and (case IV) a slot intersecting the ar-axis 
and a panel intersecting the y-axis. These conditions are 
foimd to occur when the tunnel has 4y slots with p equal to 
i (case I) or 1 (case II) or 2(2y-j-l) dots i^th p equal to i 
(case III ) or 1 (case TV). 'Since there are four sets of con- 
stants and fom different symmetrical slot configurations, it 
is to be expected that each set of constants can be associated 
with one of the symmetrical slot configurations. When 


equations (31) are set up for each of the symmetrical slot 
configurations, the following associations are found to exist: 

Case I (panels intersecting both axes): the ajj+i set has 
values different from zero. 

Case n (slots intersecting both axes): the p^+t set has 
values different from zero. 

Case m (panel intersecting the a:-axis and slot intersecting 
the 2 /-axis): the set has values different from zero. 

Case TV (slot intersecting the x- axis and panel intersecting 
the y-axis): the P^ set has values different from zero. 

Equations (30) and (31) can now be written by using only 
the constants which may have values other than zero. How- 
ever, equations (30) and (31) can be written in a more 
symmetrical form if the solution is carried a step further. 
The equations for the set of constants that are used for a 
given configuration may be written, by using tbe set ajf+i ns 
an example. 


— 1 — aiToi-hasTos'l" • • • ajj+iTo.aj+i') 


0 — 0£iTii-i-ci£3Ti8+ . . . r 


0 — cciyti+<X3yj3+ . . . ai/+i7j,sy+i J 


(32) 


where the y’s are the corresponding functions of 0, and s. 
Now consider as a matrix the coefficients of the right-hand 
side of the equations whose constant terms are equal to zero, 
and let Ai, A^, As, ... be equal to the determinants, re- 
spectively, which remain after the column that corresponds 
to the number 2j-l-l of the constant is removed. Then, 
using Cramer^s rule, 


-A, ' 

^ Acfoi — ^A.37o 3+ • • • ( — l)^Asy+i7o, ^ 

^;8 

Ai7oi — ^A.87o3-|- . . . ( — lyAij+iYo.n+i, 


(33) 


with corresponding equations for the remaining values of 


When the solution for each constant (eq. (33)) is substi- 
tuted into equation (30), it can be written 


Oj/2 


dw 


—pbT(l—s*) n V8^'~2 s* cos 20,-fl 
m (34) 


v(2*-«*)a-3*8*) W cos 20,-f-l 
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where /(z) is determined by one of the following equations, 
the choice depending upon the symmetry of the desired slot 
configuration: 

For case I, or panels intersecting both axes, 


m 


mH H—n 

J-0 

1=0 


(35a) 


For case II, or slots intersecting both axes (note that in this 
case as well as in case IV, B has the same relation to /3 as ^4. 
has to a in the demonstrated case), 

tn/4 

1=0 


For case IIT, or a panel intersecting the a;-axis and a slot 
intersecting the y-aris. 


/(2)=g (35c) 

/«o 


For case IV, or a slot intersecting the i-axis and a panel 
intersecting the y-axis, 


m+2 

m=-^ (36d) 

1=0 

The constants Bif+i, Aij, and B^j are determined from 
the foUo\ving matrices in the manner discussed in the mate- 
rial following equation (32). The matrix for case I is 


M0) 

cos 6 dd 

m 


Ch cos d dd r«i cos 3 0 dd 
J-h ji{6) J-Oi fi 


Me) 


j: 


CO»(f+l) 


dde 


Me) 


where 


mf2 


cos 20-fsO n Vcos* 0— cos* dg 

(7=1 


The matrix for case II is 

shxede f«i sin 36 d0 f®* ^ ^ ^ 




Me) Jh Me) 

hainddS 

Me) 


Me) 


(36b) 


sm 6 dB 

Me) 


where 


ni/3 


Me)=0-~^ cos 25-fs^ n Vcos* 0 ,— cos^ 6 

(7=1 


The matrix for case III is 
■ 26de 


rei ^ 

J-h Me) J-h fi' 


Me) 


p. cos (^^ 1)0 

J-^ Me) 


^ de 

h Me) 


(36c) 


dd 

, W) 


L. s 

TThere 


nf2 


/i(0)=(l— 2s* cos 204 -s^ n ■Jcoa^e—coa% 

ff-i 


The matrix for case TV is 

sin 26 dd f * sin 46 dd 


h Me) Jh Me) 

sin 20 dg 

h Me) 


^smQ+iy 

Me) 


de 



cos 6 d6 d6 

Me) 


(36a) 



sin 20 d6 

Me) 


(36d) 
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where 


ra/2 

= COS 20+s‘) n Vcos*0,— COS*0 

ff -1 


Solution for symmetrically loaded wings. — Symmetrically 
loaded wings can he assumed to be made up of vortex pairs 
of circulation —dr and dr. Since the circulation can be 
expressed as ro/(a:) where ro is the circulation at the point 
x=0, the strength of each pair becomes 

dr=ra/'(x)da: (37) 

This value of the circulation may be substituted into equa- 
tion (34) to give the contribution of each elemental vortex 
to the total flow. It is then necessary to integrate the equation 
over the entire span to obtain the flow, or 

where (dwldz)^ is the complex velocity of the nonuniformly 
loaded wing, and (dwldz)v is the complex velocity (eq. (34)) 
of the uniformly loaded wing with x substituted for s and 
a circulation of imity. 

Corrections for interference due to lift. — ^The interference 
complex-velocity fleld at the far-downstream position is 
determined by subtracting the complex velocity of the free 
field from the complex velocity of the constrained field, or 


obtain the correct drag is 

AO^=S^ (42) 

In the case of the slotted tunnel, it is convenient to use the 
quality factor k and then express the two corrections as 

(43) 

and 

AOn=kb ^ (44) 

where 5 can be determined from the hterature on lift inter- 
ference of closed tunnels. 

Solution and interference quality factors for tunnels with 
cross sections other than circular. — ^The solution of the wall- 
interference problem for tunnels having cross sections other 
than circular can be obtained by following the method of 
paragraph 14.6, reference 8. This method requires that a 
function 2 =/(f) be found which will conformally transform 
the interior of the tunnel cross section in the 2 -plane into 
the interior of the tmit circle lfl=l in the f-plane. It is 
also necessary thaty'(f) does not vanish or become infinite 
in the unit circle |f|=l. The complex velocity of the 
tunnel in the 2 -plane is then 


u~w= 


dw dt 
df dz 


(46) 


dwi_dw ■ isT 
dz dz 


(39) 


A useful parameter which indicates the effect of slotting a 
tunnel is the ratio of the interference of that tunnel to the 
interference of a closed tunnel. This parameter, which is 
the same at the far-downstream position as at the model 
and which will be called the quality factor k, can be expressed 
mathematically as follows: 


k= 


dz 

ir(z*— «*) 

dwc 

, isF 


(40) 


where dwcldz is the coipplex velocity of the closed- 
tiumel configuration. 

Once the values of k are determined as a function of the 
semispan s for a specific tunnel, the interference can be cal- 
culated. It is shown in reference 7 that the interference on 
the lift can be expressed as an angle which may be added 
to the measmed angle of attack to obtain the true or free- 
flight angle of attack. Reference 7 gives the correction 
angle in radians as 



(41) 


where 5 is a factor which is determined from the geometry 
of the tunnel configuration, 8 is the wing area, O the tunnel 
area, and Oi, the lift coefficient. Also from reference 7, the 
increment which must be added to the measured drag to 


In the function dwjdt of equation (46), the B/s which de- 
termine the slot edges of the f-plane are the transformed 
values of the slot edges in the s-plane and the values of i! in 
the f-plane are also determined from the transformations of 
the points s in the 2 -plane. Once these values are deter- 
mined for the transformed tunnel in the f-plane, the velocity 
field dnvjd^ in the f-plane may be computed. The velocities 
in the z-plane may then be computed by equation (46). 

The interference and the quality factor may now bo de- 
duced from equation (46). Subtracting the free field from 
the complex velocity given in equation (46) gives for the 
interference velocity 

dwi dw df . isF / . 

dz “df d2'^v(z“-0 ^ 


The interference for any tunnel can be determined from this 
equation; however, in many ca^es, the use of a qualify factor 
may be more convenient. The quality factor for this class 
of tunnels can be written 


■isF 


df dz 

v(2*— s’) 

dwo df , 

isT 


(47) 


where dwo/df represents the complex velocity of the trans- 
formation of the closed tunnel in the f-plane. 

A simplification of equation (47) can be made in case the 
transformation y(f) may be approximated by cf for points 
near the center of the tunnel. If s is sufficiently small, its 
value in the f-plane may therefore be represented by 8— cl, 
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where I represents the distance to the point which locates 
the transformation of the tip vortex in the f-plane. Sub- 
stituting these approxiniations into equation (47) gives for 
the quality factor: 

dw 1 . idV dw . UT 

. dr ..ox 

dWa 1 , iclT dWc , fir ^ 

Equation (48) shows that the quality factor near the center 
of the original tunnel in the z-plane is approxunately equal 
to the quality factor of the corresponding circular tunnel 
in the f-plane. Thus, to obtain the interference, only the 
quality factor of the circular tmmel needs to be computed, 
provided the interference of the original closed tunnel is 
known and the approximation z=cZ is valid throughout the 
region in which the model is located. 


APMJCATIONS OF THEOEY TO VAEIOUS TDNNEin 

Two slots located symmetrically across the x-axis. — ^Tun- 
nels which have two slots located symmetrically across the 
2 - axis (fig. 2(a)) are treated under case IV. Thus, their 
complex velocity field may be expressed by equation (34) 
and equation (35d) where y takes the values 0 and 1. Since 
the integrals across the two panels are identically zero, the 
matrix of the integrals has no meaning and hence is not used. 
With these considerations the complex velocity may now 
be written 


dw 

dz 



The wall interference is determined by subtracting the 
complex velocity for the free field from equation (49), or 

dwi_ isV r (1— 2 ^)Vs*— cos 2^1-f-l 
dz~ 7r(z“— s*) L(l— z"s^)V 2 *— 2z* cos 201 - 1-1 J 

The quality factor k is determined by dividing equation (50) 
by the closed-tunnel interference which is 

^ isv ra-«‘)(i+25„ii 


r (i— s^) 

ir(z*— S*) L 1~ 


(1— 2a^ cos 2gi-|-l , 

, (1— s*s*)Vz^— 2z^ cos 201 - 1-1 , 

This function can be written 

cos 201 - 1 - 1 — (l—g^8^Vg^—2g^ cos 2gi-fl 
(g“— s*) -y/?— 2?^os"25rTd 


The effects of the slots on the interference along the span of 
the model are obtained by substituting x for g; then, 

(1—2^ 78^—28^ COB 2gi~|-l— (1— x^a^-y/x^— 2x^ cos 2gi-pl 
( 2 ®— 8 ^-y/ 2 ^~ 22 ^ cos 25i-l-l 




(a) Case IV; 

2 slots, p=l. 

(c) Case I; 

4j slots, p=i. 

(e) Case III; 

2(2j+l) slots, p=t. 
(g) Single slot. 


(b) Case III; 

2 slots, p=t. 

(d) CaseU; 

4j slots, p=l. 

(f) Case IV; 

2(2j+l) slots, p=l. 
(h) Single slot. 


Fiqurb 2. — Various symmetrical slot configurations investigated for the 
lift interference. 


If 2 and 8 are sufficiently small][so that the approximation 


2 * . o , 


/2^— 22® COS 201-1-1=1—2® COS 20i-f~ sin® 26i 


is vahd, then 

, 2 cos 2 ^ 1 — 22 ®— ( 2 ®-f «^'*3in® 26i 

2 — 22 ® cos 20 id- 2 * sin® 26i 


(55) 


The quality factor is seen to be insensitive to values of either 
2 or 8, so that for this tunnel with relatively small models 
(s=0.25 or less) the effect of slotting the tunnel on the lift 
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of any model may be obtained by multiplying the closed- 
tunnel interference for that model by the slotted-tunnel 
quality factor. 

Pistolesi (ref. 5) shows that the interference quality factor 
for this tunnel with a vortex doublet in the center is, in the 
notation of this report, cos 20i. Equation (65) reduces to 
the same value when x and h are set equal to zero. 

Two slots located symmetrically across the y-axis. — ^Tun- 
nels which have two slots located symmetrically across the 
y-axis (fig. 2(b)) are treated under case IH (see eq. (35c)). 
Thus, their complex velocity may be expressed by equation 
(34) with/( 2 ) given by expression (35c), so that it can he 
written 


dv) — isr(l— 2s^ cos 26 i+ 1[2 Aqi^—A 2(1 +z^] 

dz — s*)(l — — 22 * cos 2 fli-t-i[2AoS^— .Aj(l+s0] 

(56) 


where and A2 are determined from the matrix (36c) and 
are 



cos 28 dB 

(1 — 2s® cos 26-{-s*)-yJcos‘d—cos^6i 

de 

(1 — 2s® cos 20-f-sO V cos®0— cos®fli 


(67) 

(58) 


If 2 and s are considered small, the complex velocity of 
this field may be approximated by 

dwj •isr (AjS®— 2Aq)-h(l— s^(j4.a— 2.4 o 2®) cos 201 , . 

~d^~ (^a-2^oS*)(l-2®s®) ^ 

from which the quality factor becomes approximately 

, (AaS®— 2A(Q-K1— S^(Aa— 2 iAo 2 ®) cos 20i 
~ A2—2AoS^ ^ ^ 

which reduces for z=s=0 to 

OA . 

*=-^+cos20i (Bl> 

This quality factor is the same as the one presented in 
reference 5, with the necessary change of sign to adapt it to 
the notation of this report. 

Tunnel with ij slots located symmetrically with respect 
to the X- and y-axes. — 'T unn els with slots located sym- 
metrically with respect to the x- and y-axes can have either 
the symmetry of case I (eq. (35a), fig. 2(c)), if panels inter- 
sect both axes, or of case II (eq. (36b), fig. 2(d)), if slots 
intersect both axes. 'Thus, the complex velocity can be 
expressed by equation (34) with/( 2 ) given by either equation 
(35a) or (35b), depending on the desired slot configmation. 
As an example, the complex velocity for a tmmel containing 
four slots with panels intersecting both axes may be written 


—i8T(l—8*)(jl^8*—28‘ cos 20,-fl^[Ai2*(2®-l-l)— . 43 ( 3 ®-+-!)] 
v( 2 ®-s®)a - cos 20,-f-l)[4i8®(s®-l-l)-^,(8'’-f 1)] 


(62) 


whore Ai and A3 are determined from equation (36a) and 
are expressed as 


A,-p- 

J-hn- 


cos 30 d0 


(1— 2s® cos 20-t-s*) n Vcos®0— cos*0 


and 


-®i n — 


cosdde 


(1 — 2s® cos 20-fsO n VcO3®0— COS®0, 


(63a) 


(63b) 


If the tunnel contains eight slots with panels intersecting 
both axes, 


dw 


i3T(l—8*) nV«^— 2s* cos20,-l-l 
0=1 X 


4 - . - 

7T ( 2 ®— s®)(l — 2 ®8®) n 22® cos 20j-}-l 

B-1 

Ai2*(2®+l)-.A32®(2«-H)-hA(2^°+l) 

AiS-‘(s®-M)-A[3S*(s“+1)+A(s'“+1) 


(64) 


where Ai, .4s, and A^ are equal to the determinants remaining 
when the first, second, and third columns, respectively, are 
removed from the following matrix: 


r®i cos odd 

cos 30 dB 

r«i cos 58 dB~ 

fM 

J-h Me) 

J-oi M^ 

r®3 cos 0 dB 

pj cos 30 d0 

cos 50 dB 

r~ 

jh Me) 

Jh Me) J 


(65a) 


where 


/i(0)=(l— 2s® cos 20-}-8*) n Vco3* 0— COB® 0, (66b) 

0-1 

If the tunnel contains 12 slots with panels intersecting 
both axes. 


, — isr(l— 8^ n Vs‘— 2s® COB 20, -fl 

dw g=i 

dz 


-X 


7T ( 2 ®— 8®) (1 — 2 ® 8 ^ n 2z* cos 20(,-l-l 

0-1 


4i2«(2®-H)-4L3g^(z»-H)-h4»2®(4°-j-l)-47(2»^-H) . . 

AiS«(8®-|-l)-488‘(8»+l)-l-43S®(s‘''-l-l)-4r(s>^+l) ^ ’ 


where Ai, A3, A3, and A^ are equal to the determinants 
remaining when the first, second, third, or fomth columns, 
respectively, are removed from the following matrix: 


pi cos odd 

r»i cos 38 dB 

r*i cos 5BdB 

pi cos70d0" 

J-h Me) 

J-h Me) 

J-e. Me) 

J-Oi Me) 

p3 COS BdB 

p3COs 38 dB 

p 3 COS 56dB 

f«3COS 78 dB 

J03 Me) 

Me) 

Me) 

Jsi Me) 

f»5COS BdB 

pscos 38 dB 

p3 cos 5Bd6 

pt COS 7BdB 

1 

Je, Me) 

J03 Me) 

Jh Me) J 


(67 a) 


where 

/i(0)=(l— 2s® cos 20+8^ n Vcos® 0— cos® dg (67b) 

0=1 
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If the slots are of equal "width and evenly spaced, equations 
(62), ("64), and (66) can be simplified by substituting the 
terms 

VI-22" cos 201+2*" 
and 

Vl— 2s®* cos 25i+s^ 

for the multiple products in equations (62), (64), and (66), 
and ^ 

V 2 n 9 A 
cos* y 0— cos* Y 

for the multiple products in equations (63a), (63b), (65b), 
and (67b). With these substitutions the limits of the 
integrals in the second row of the matrix (66a) are changed to 

^—01 for the lower limit and for the upper limit, and 

those of the second and third rows in the matrix (67a) are 

changed to ^ — 0i and §—6i for the lower limits and ?+0i and 

^+01 for the upper limits. 


The interference in these tunnels is obtained by computing 
the velocity and subtracting the free-field velocities, then 
dividing by )( since the computed fields are for a great 
distance downstream and (see ref. 6) the interference veloci- 
ties at the model are one-half those downstream. 

Tunnels with 2(2y-|-l) slots located symmetrically with 
respect to both axes. — ^Tunnels -with 2(2y-M) slots located 
sjrmmetrically with respect to the x- and y-axes can have the 
symmetry of slot location of either case III (expression (35c), 
fig. 2 (e)), with a panel intercepting the x-axis and a slot 
intercepting the y-axds, or of case IV (expression (35d), 
fig. 2 (f)), with a slot intercepting the x-axis and a panel 
intercepting the y-axis. Thus, the complex velocity may 
be expressed by equation (34) "with f(z) given by either 
equation (35c) or (35d), depending upon the desired con- 
figuration. As examples of this class and of the simplifi- 
cation due to evenly spaced slots of equal "width, the complex- 
velocity fields for tunnels contaioing six evenly spaced 
slots of equal "width can be "written for case m, where 
panels are on the x-axis. 


dw -isr(l-a^Vsi*-28»cos60i+ l2AoZ*~A^i^(z*+l)+Ai(2^+l )] 
dz t( 2 *-s*)( 1 - 2 * 8 *) V 2 “- 22 « cos 601-M [2AoS<-AjS*(8^-M)+A,(s«-i-l)] 


where A^, Aj, and At are evaluated from the following matrix 
by the same methods used in equation (64): 


r«i de 

'**1 COS 20 dd 

'**1 cos 40 dd 


Me) J 

-«i Me) 


cos 20 dd 

'’'7"’^ cos 40 dd 

3 " 

i-e, Me) J 

«|*t 

3* 

1 


(69a) 


where 


/i(®)=(l— 2«* cos 20-fs^ Vcos* 30— cos* 30i (69b) 

For case IV (fig. 2 (f)), where slots are on the x-axis. 


d-w -isr(l~s*) V8>*-28'» cos 601-H 

dz ^(2*-S*)(l-2"s*) Vz‘*- 22« cos 601+1 [5aS*(8^-l)-.B4(8«-l)] 


where and are determined from equation 36 (d) and 
are expressed as 


and 



sin 20d0 

(1—28* cos 20-1-8^ Vcos* 301— cos* 30 


(71a) 



sin 40 d0 

(1—28* cos 20-f-s^ Vcos* 301— cos* 30 


(71b) 


The interferences are determined as m the pre"rious 
sections. 

Equations (68) and (70) may be extended to 10, 14, . . . 
4j-]-2 slots in the same manner that equation (62) was ex- 
tended to include 8 and 12 slots. 

Single slot located symmetrically across the x-axis. — Since 
tunnels which have a single slot located symmetrically 
across the x-axis (fig. 2 (g)) contain an odd number of slots, 
summation (17) rather than (11) must be used in equation 
(14), which, when expressed for a single slot, becomes 


dw 

dz 


■■±.0 
n«0 \ 


On COS ^ ^+ln S1H 


) 


(s* S*)(l 2^8 * ^QQS d— cos 01 


Equation (72) may be re"written, "with several modifications 
after substituting z=e^, as 

dM?_goz(l-|-g)-|-(ii(l-|-g^-f-i[6og(l — g)-f-6i(l — g^] 

dz V2 (z*-8*)a-2*8*) Vi='2g cos 01-1-g^ ^ ^ 

The constants in this case must be evaluated from the 
circulation about the poles at s and — s and the nonexistence 
of sources at those points, as the potential condition is 
automatically satisfied because the potential is constant 
across a single slot. The circulation and noneSdstence of 
sources (eqs. (27) and (28)) may be expressed about the 
pole 8 as 

j. [ 27ri{fflos(14-£)+gi(l+8^+t[6os(l— 8)-b6i(l— 8*)]} 

V2 2s(l — 8^ Vl— 28 cos 01-1-8* 

and about the pole — s as 

P I 27rz{ao( 8)(1 — 8)-f-ai(l-|-8^-|-f[6o(— 8)(l-|-8)-h6i(l — s^ ] } 

V2(-28)GL-8Vi+2s cos 8i-1-s^ 


2 


(72) 


(74b) 
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If there are to be no sources at the poles, a- must equal zero. It may be seen from inspection that Oo and Oi are equal to 
zero and that bo and bi have solutions other than zero. When these solutions for bo and bi are substituted into equation 
(73), it becomes 

dw —isT ('r , ^ ^ , , 

Vl-2z cos \ (1+s) V1-2S cos (1-s) ^l+2s cos 0,+s^J e(l-z)+ 

[— (1— s®) Vl+2s cos 0i+s^+(l+s®) Vl~2s cos &i+s*J (1— ^ (76) 

The interference velocity is determined as before by subtracting the free-field velocity hence, the inter- 

ference complex velocity is 

dwr —isT C r , : , ^ , 

(1-e^s*) Vl-2z cos 2(1-^)+ 

[( 1 -fs^ Vl~ 2 s cos 0i4-s^— ( 1 — Vl+ 2 s cos 5i4-s^] ( 1 — s®)— 2 s ( 1 — 2 ®s*) ■^jl—2z cos fli-f (76) 

The quality factor k is determined by dividing equation (76) by the dosed-tunnel interference velocity, which is 


— isr/ir(l— 2 ®s®). After performing this operation and combining 


r\-f +V*A 


^‘~ 2 s (g®— a®) V l— 2 z cos [(!+«) 0-—^ ( 2 +«) (X+zs) ■y/l—2s cos 0 i+s®— ( 1 — s) ( 1 — s) (z—s) (l~zs) Vl+ 2 s cos Oi+s^— 

2s (1 — 2 ®s*) Vl — 2z cos 01+2®] (77) 

This factor can he checked by letting 0i=O°, which represents a closed tunnel, for which k can be shown to be equal to + 1 ; 
and if 0i=18O°, which represents the open tunnel, then k can be shown to be equal to —1, 


Equation (77) can be considerably simplified by using the 
following approximations, in which z (or s) is assumed to be 
small with respect to 1 : 

i/l—2z cos 01+2^=1 — Z cos 01 + 2 (78) 

and 

Vl+ 2 s cos 0 x-fs*=l-f s cos 01 +s® ( 79 ) 


After multiplication and collection of terms, equation (77) 
simplifies with the use of equations (78) and (79) to 


. sin* 01 
cos 01 = z 


k=- 




1 — Z cos 01 + 3 * 


sm* 0 i 


(80) 


If terms of the order s* and higher are eliminated, equation 
(80) may be approximated with 

k = — i { 1— 2 cos 01 — cos* 01+3 [ 2 s*+( 2 — cos 00 sin* 0 j } 

(81) 

Equation (81) shows that the interference of this tunnel has 
an odd function component along the a:-axis which will 
produce a variation in the effective angle of attack that will 
cause the model to have a rolling moment. Since extension 
of flow fields of this type (e^ odd number of slots symmetri- 
cally located with respect to the a;-axis) to greater numbers of 
slots win not eliminate the odd powers of 3, the rolling 


moments will continue to exist. It is therefore to be con- 
cluded that tunnels containing an odd number of slots 
symmetrically located with respect to the a:-axis introduce 
an extraneous moment into the data. 

Single slot located symmetrically across the y-axis. — For 
tunnels which have a single slot located symmetrically 
across the y-axis (fig. 2(h)) summation (17) is again used 
rather than ( 11 ) becaiise of the odd number of slots. This 
summation when used in equation (16) results in the proper 
slot symmetry, and the complex velocity can therefore bo 
expressed as 


dw z aocos-+aicos^+6osm2+6ism-^ 

d3~(3*— s^(l— 3*s*) Vsin t?— sin 0i 

(82) 

Since z—e^, the complex velocity may be rewritten after 
substitution and combining various terms as 

dw^ (X + i) [aog(l + g) +ai ( 1 + 3*)+ 6 oi 3 (l —z)+ 6 ri(l — 3*)] .gg. 

dz 2(3*— «*)(!— 3 * 8 *)Vl~ 2 i 3 sin 0 i— 3* 

The constants Oo, ai, bo, and 61 must now be evaluated 
from the circulation about each of the poles and the condition 
of continuity. The two conditions about each pole may bo 
expressed as 

^ 27^^(l+^)[go8(l+8)+gl(l+8®)+6o^a(l— a)+&ii(l— a*)] 

4s (1 — s^ Vl~2'is sin 0i— s* 

(84a) 
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about the pole s and 

p I ^•^_ 27ri(l+i)[— gos(l— s)+ai(l— 6oZ8(l+8)+5i^(l+g^] 

—4s (1— s^Vl+2w sin 0i— s* 

(84b) 

about the pole — s. The strength tr must equal zero in order 
to satisfy the condition of nonexistence of sources. In 
order to simplify the rationalization of equations (84), the 
following substitutions are made: 


C7+iZ?=Vl— s^+2ts sin 6i 


Also let 


0—iD=-\jl—8^—2is sin 01 

yli=aoS (1 +s) +Oi (1 +s^ 

.d.2= — OoS (1 — s) +<ii (1 — s®) 
5i= 6oS (1 — s) + 6i (1 — s®) 

Bi— — 6os(l+s)+6i(l+s®) 


(85a) 


A=- 


2sr(l— sQVl— 2s® cos 2di+8* 

TT 


(85b) 


(86) 


The four simultaneous equations (84) can no"w be written, 
after rationalizing the denominator, 

4=R.P. (1-tO {G+iB) (ALi+iBi) 

0=I.P. (1-i) {G+W) (^i+iBi) 

A=R.P. (1-tl (C7-iZ))(A2+iB2) 

0=1. P. (1 {G~W) (udj+iBj) 

which may be restated as 

A=A,{G+B)+B,{G-B) 

Q=-A,iG-D)+B,{G+D) 

A=Ai(P-B)+B^{G+B) 

0 = -A3 (G+B) +Bs(G-B) 


(87) 


If Ai and Bi are determined from the first two and A3 and 
B3 from the last two of equations (87), it may be seen that 


'^1=<*0S (1 +s) +Oi (1 +8^ 
-OoS (1 — s) +Oi (1 -s®) 


(88a) 


^ jS^p^3) =-8i=^os(l-s)+f.i(l-s®) 

2 ~ 60S (1 +s)+ &i(l +s®) J 


(88b) 


Equations (88a) "when solved for Oo and Oi and equations 
(88b) -when solved for 60 and bi give 


A (G+B)(1-8^-(,G-B)(1+8^ 


^~<7®+Z>® 

4s(l— s^ 

A 

(a+Z?)s(l-s)+((7-X>)s(l+8) 

^ 0®+Z>® 

4s(l— sO 

A 

(G+B) (l-^-(G-B) (1+s®) 

6 ’®+Z>® 4s (1-8^ 

A 

((7+Z))8(1-8)+((7-Z>)s(1+s) . 

" (7®+Z>® 

4s(l-s«) ^ 


Equations (89) may be simplified to 


Uo — — ba — 


2A(B-G8<^ 


cti — bi — 


4(C®+I?^s(l-s^ 

2A8 (G-Bs) 

4(C“+i)®)s(l-s^ 


(89a) 
(89b) 
cEo (89c) 
(89d) 

(90a) 

(90b) 


Letting ao=— &o and. (h=bi, the complex velocity equation 
(83) may be rewritten 


dw 


2[g(ao+aig®)+2(gos®+ai)] 


dz 2(2®— 8^(1— g®s^Vl~2fg sin 2® 

Substituting for Oo and Oi from equations (90) gives 


(91) 


d'u;_ 2A { 2 [P — (7s®+s {G—B8) 2 ®] [s (O—Bs) +(B— (7s®) s®] } 
dz 4((7®+Z>®)s(l— 8^(2®— s®)(l— 2®s®)Vl— 2i2 sin Oj— 2® 

(92) 

Now it may be shown that 


A=-— (l-s«)((7®+Z>®) 


(93) 


so that 

dm 
dz 


r 2[(P-C's®)+s(g-Ps) 2®]+i[s((7-J>s)+(J>-gs®) 2®] 

7T (2®— 6*)(1 — 2®6®)V1— 2-12 sinOi— 2® 

(94) 


where 


(7+iIl=Vl"“S®+2fs sin di 


The interference velocity is obtained by subtracting —is r/ir(2®— s®), so that 

dwi 


isV C 8jG-B8)+(P-G8^^-iz\B-G8^+8iG-B8)^ 
dz tt( 2® s®)"^ g(l — 2*8®) Vl — 2i2sin0i— 2® 


} 


( 95 ) 
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Now divide by the dosed-tunnel interference to obtain the quality factor, so that 


or 


—far a((7— f7s^+s(<7— Z?s)g^— s(l— g^a^yi— 2i2smgi— 

7t( 3^— a(l— 2^g*)Vl— 2i2Bmgi— 

— isr r (i— 8^(i+g^ .~1 

ir(2*— L 1 — 2*8* J 

8(g-Z>8)+(P-Cb*)2*— i2 [Z>-Cfe^+8(g-Z>8)2*]-8(l-2*8*)-^^-2^2Singl-i^ 

8 ( 2 * — 8 ^Vl — 2iz sin di — 2 * 


( 96 ) 


(97) 


The interference quality factor at the center of the tunnel 
may be examined by letting 2 = 0 ; then 


b 


G—Ds — 1 
-8^ ■ 


(98) 


From the evaluation of the constants, equation (94), 

G=Xj\ — 28* cos 201 + 8 ^ cos ^ (99) 


Once the constants, K, K', and m are determined, the slot 
location and the half-span length I may be computed. 
First, consider slots which are located on the top and bottom 
of the tunnel (see fig. 1(b)). In reference 8, paragraph 14.8, 
it is shown that the top of the tunnel may be expressed by 
ih 

z=x+-^ or \z=\x+iK'. It is also shown that the trans- 
formation (106) may also be expressed 


where 


Z?=Vl — 28*cos 201 -f 8^ sin ^ 


tan-i 


28sinfli 

~T=^ 


aoo) 

( 101 ) 


When 8 is small, equation (98) may be simplified by first 
aI>proximating O and D. These approximations are, for 
small values of 8, 

(7=1 (cos 25i-l-sm* 0i) (102) 

D=s sin 0-03) 

When the values for O and D given ia equations (102) and 
(103) are substituted into equation (98), it becomes 

, ... cos * 01 

K=sm0i-j — 2 — (19^) 


This equation also checks with the residts given in reference 
5, with the usual change in sign to conform with the notation 
of this report. 

lift interference of rectangular tunnels. — The function 
which transforms a rectangular tunnel into a circular timnel 
is given in paragraph 14.8, reference 8, as 


^= 


X 2 , X 2 
cn^ 


(105) 


, \z j \z , 

where sn » dn and cn 

of 2 and X/2 is defined by 


X 2 

~2 


are Jacobian elliptic functions 


\_K_ K' 
2 a h 


(106) 


where K and 'iK' are the quarter periods and a and h are the 
breadth and height of the tunnel. (See paragraphs 14.7 and 
14.8 of ref. 8 for further information concerning these func- 
tions.) 

If a and h are given, then K, K', and m, the squared 
modulus (see ref. 8, paragraph 14.8), are uniquely determined. 




1— cnX 2 
l-fcnX 2 


(107) 


so that the equation for the top of the tunnel becomes, in 
the f-plane. 


^ 1 — cn(Xr-fi7irO 
1 -j- cn(Xa:-t-iK'0 


(108) 


From reference 8, paragraph 14.8, cn(Xa-l-iiir') is equal to 
— m‘^dsfXa;, which wiU, for convenience, be called in. Now, 
on the circle. 


f*=8*«= 


1 — in 
l+i^ 


(109) 


and if the real and imaginary portions are separated, 

1— /J* m— (dsXa;)* 


cos 20= 


or 


sin 20= 


1-1- /n* m-)-(d8Xa:)* 

2n 2dsXr 1 

1+m"~ Vm 2 I (dsXr)' 


( 110 ) 


( 111 ) 


m 


The 0/8 of the circular tunnel in the f-plane are determined 
from the above equations by the locations of the slot edges 
x, in the 2 -plane (see fig. 1(b)). 

If the slots are on the side of the timnel, a similar analysis 

starting with z—^+iy or 'Kz=K+i\y wfil show that 


• nff 2iVl— w sd (iXy,m) 2^jl—m sd (\y, 1— m) 
sm (^Xy^^)]2-l+(l_m) [sd(\y, 1-m)]* 

( 112 ) 


As before, equation (112) may be used to determine the 
0/8 when the slot edges are functions of y. 

The value of ( is given in reference 8 as 


„ 1— cnX6 

l-l-cnX6 


(113) 


With I and the d/s known, the complex velocity of the tunnel 
in the f-plane may now be computed, and the interference 
may be determined from equation (46). 
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The calculation of the interference for rectangular timnels 
may be simplified by applying the conclusion following 
equation (48), which states that the quality factor of the 
transformed tunnel is approximately equal to the quality 
factor of the original tunnel, provided the transformation 
function can be approximated with z=cf. Equation (105) 
is shown in reference 8, paragraph 14.8, to be approximately 
equal to f=\ 2 / 2 ; hence the conclusion following equation 
(48) is valid for rectangular tunnels. Thus, it is necessary 
only to determine the quality factor for the transformed 
tunnel and to apply it to the correction for the fully closed 
rectangular tunnel (see ref. 8, paragraph 14.8, for this 
correction) to obtain the correction for the corresponding 
slotted rectangular tunnel. 

Circular tunnels having general slot configurations. — 
Equation (16) is used to obtain the complex velocity for a 
configuration which contains a general slot distribution. 


Equation (16) may be vpxitten 


d/m 


“pzich cos t?+6i sin t?) S (“»+i/3n) cot"^ 


n-0 


dz 


I A ft 
(S*— S*)(l — n -YCOt^— COt^ 


(114) 


where p has the same definition as given for equation (25). 
If 2 is substituted for by using the relation z=e^, equation 
(114) becomes 




(2®— (1 — 2*S^ cot ^ 

Equation (115) may be reduced to 

p[oi(l+2*)+ iSi(l— 2^] S(a»+il3,0(2— l)"~’'(z+l)V 


(IIQ 


dz 


2m I ^ 

2 ( 2 *— 8^(1— 2 *s*) n i(z+l)—(z—l) COt^ 

0=1 V -li 


(116) 


The constants Oi, bi, oc, ai, . . . otn, and fio, Pi, ■ P„ are determined in the same manner as was used for the sym- 
metrical slot configurations; that is, the circulation about both poles is equal to r, the source strength a is zero, and the 
potential in each slot is equal to zero. The conditions on the circulation and source strength may be determined from the 
line integral about the poles s and —a. For the pole at a, 


and for the pole at —a, 


2Hp[<h(l+a^+ibiO—a^] X) (a,-l-ij3,)(l-s)’»-»a+s) V(-l)"-’* 

T-{-lcr= — 

2m I 

4s(l— n -y t(s-l-l)— («— 1) cotj 
27rip[ai(L+a^+ibi(l—3^]^(an+ip„)(l+a)’’-’'(l—8)”i"(—l)”-" 


r-|- 'io'= 


2m j 7 

4s(l— 8^ n .»/t(l— s)-f(l-hs) cot-; 
ff-1 V 


(117a) 


(117b) 


The condition on the potential in each slot may be determined from integrating the velocity over each panel, or 

d~\ 

(oi cos 6+bt sin ^2j(«i+i/5i) cot"^ J<^ 


>=f ‘ 

Jpantl 


I w iT 

(1— 2s* cos 2e+a*) n '»/ COt ;5 — C 0 t-;f 

t-1 \ ^ 


(117c) 


The constants Oi, bi, oo, ai, . . . a,, and fio, Pi, •• • Px 
are determined from the solution of the set of simultaneous 
equations (117). Once the constants are determined, the 
complex velocity may be determined, and the interference 
may then be computed in the same manner as was used for 
the various symmetrical cases. 


RESULTS AND DISCUSSION 

INTERFERENCE QUAUTY FACTORS FOR SEVERAL CIHCOLAB TUNNELS 
WITH SYMMETRICALLY LOCATED, EVENLY SPACED SLOTS 

The quality factors for circular tunnels containing various 
symmetrically located, evenly spaced slots of equal length 
and wings of very small span (s->0) are given in figure 3. 
The curves for tunnels that have panels across the z-axis and 

413072 — 57 25 


contain 4, 8, or 12 slots are calculated from the formula 

cos ^^+1^ ® ^ 

C03*^ 

— — k- V 

h COS ( "2 — ^jOdB 

2 ^ o 7^1 
cos* Y 0— cos* 

where 9 is defined in figure 1 and m is the number of slots. 
Equation (118) may be derived by considering the value of 
equation (34), expressed for symmetrically located, evenly 
spaced slots of equal width, when a and z are equal to zero. 
It is also the negative of the relation given in reference 5 for 
the same tunnel configurations, the minus sign being used to 
conform with the notation of this report. 
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Number of 
slots 


Slot 

configuration 



Percentage of tunnel wall open 


Number of. Percentoge of 
slots tunnel wall open 


Tunnel cross 
section 


Slot 

configuration 



Ratio of wina soon to tunnel dlometer. b/D 


Figubb 3. — Quality factors for small-span wings (s-*0) in circular 
tunnels containing various numbers of evenly spaced slots of equal 
length. 


Figtob 4. — Variation of the quality factor at the center of the model 
with the span of the modeL 


An analysis of figure 3 shows that, for all configurations 
except the one with two slots across the x-axis, only a small 
percentage of the tmmel wall must he opened in order to 
obtain no interference at the center of a small model and 
that the amount of opening required rapidly becomes smaller 
when the number of slots is increased. It is also noted that 
the change in quahty factor is very rapid at the nuU-iater- 
ference condition, so that an accurate estimate of the iuter- 
ference is difficult. 

The importance of slot location in determining the quality 
factor is indicated by the large divergence of the two curves 
for tunnels containing two slots. It can be shown, however, 
by comparing the quality-factor functions for the conditions 
of slots across the x- and y-axes and panels across the x- and 
y-axes, that as the munber of slots increases, the quality 
factors wiU approach each other. Since the quahty factors 
approach each other, it may be expected that the quality- 
factor curves for tunnels containing 8 or 12 evenly spaced 
slots of equal length will be approximately correct for any 
slot location. Thus, the quality factor should be approxi- 
mately correct, even though the model is rolled in the tunnel. 

VABIATION OF QUALITY FACTOR WITH SPAN OF MODEL 

The effect of model span on the quahty factor k at the 
center of several slotted tunnels is shown in figure 4. Exam- 
ination of the curves shows that the quahty factor is fahiy 


constant over an appreciable range of spans, in that it does 
not vary more than 10 percent for spans of 0.5 to 0.6 of the 
tunnel diameter or width. The quahty factor for the tmmel 
containing two slots is seen to vary more, in that the quahty 
factor changes 0.10 for a change in the span from zero to 0,46 
of the tunnel diameter. The quahty factor of the single slot, 
however, is seen to remain fairly constant regai'dless of span. 

The spanwise variations of the quahty factor for a model 
with a span of 0.75 of the tunnel diameter is shown in figme 5. 
This figure shows that the quahty factor for the tunnels con- 
taining 8 or 12 slots is about 0.3 larger at the wing tips than 
it is at the center, whereas for the tunnel containing two slots 
it decreases about 0.6. These changes indicate that the 
spanwise change in k can depend upon the bovmdary condi- 
tion at the intersection of the waU and x-axis; that is, if a 
panel intersects the x-axis, k wfil become larger, and if a slot 
intersects the x-axis, k wfil become smaller. Such a tendency 
would be difficult to prove, however, so that an analysis of 
the variation in k for any particular tunnel would require a 
computation of k along the span for that model. These 
variations in k also indicate that the lift corrections for a 
model with a span as great as 0.75 of the tunnel diameter can 
be roughly approximated by using the value of k at the center 
of the tunnel; however, a more accurate correction would 
involve the use of an average value for k as well as an average 
of the load over the span. 
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Number of Percentage of Slot 

slots tunnel wall open configuration 



Ratio of distance along model semispan 
to tunnel radius, x/R 


Fiqubb 6. — Variation of the quality factor along the span in several 
tunnels containing equally spaced, symmetrically located slots of 
equal length. 6=0.76D. 


A r6suia6 of the observations made from figures 4 and 5 
indicates that, for these specific slot configurations at least, 
the quahty factors for tunnels containing 8 or 12 slots, pro- 
vided they are symmetrically located, equally spaced, and 
of equal length, are more nearly constant throughout a 
greater portion of the central region of the tunnel than are 
the quahty factors in tunnels containing two slots. The 
fact that the quahty factor is moderately constant even at 
spans of 0.76 permits the correction for lift interference to 
be made by computing the interference of the closed tunnel 
having the same cross section and then multiplying that 
interference by the quahty factor for the center of the slotted 
tunnel. 


QUAUTY FACTORS m TUNNELS CONTADUNG A SmGLE SLOT 

The variation of the quahty factor with percentage of 
waU open in a tunnel containing a single slot symmetrically 
located with respect to the a;-axis is shown in equation (80) 
and in figure 6 to be somewhat different from those previ- 
ously studied in that k is greater on the half-span of the 
^ving which points toward the panel and is smaher on the 
side which points toward the slot. This variation of k means 
that a spanwise variation in angle of attack exists which 
would cause the model to roU. Thus, a tunnel of this type, 
that is, a single slot symmetricahy located with respect to 


Point at which interferencb Slot 

Is calculated, x/R configuration 



Percentage of tunnel well open 


FiauBE 6. — Quality factors for a tunnel containing a single slot tvith a 
model of span equal to one-half the tunnel diameter. 

the a:-axis, has interferences which are more diffi cult to cor- 
rect for because of the introduction of an xmnecessary roU- 
ing moment into the measmed data. 

This rolling moment does not exist, however, if the slot 
is symmetrically located with respect to the T/-axis. It is 
shown in equation (97), however, that the induced velocities 
wiU not be normal to the span but will have a small compo- 
nent along the span. Such a component should not affect 
the total forces seriously, since it would not affect the total 
induced velocity. However, the effects of this component 
should perhaps be considered in the treatment of load dis- 
tributions along the span. Thus, the pecularities of the in- 
terference effects of timnels containing a single slot symmet- 
rically located with respect to either axis indicate that its 
effects will be more difficult to correct. 

If the equations for the single-slot case are extended to a 
larger odd number of evenly spaced, equal-width slots, it 
can be expected that both the rolling moment and the cross 
flows will become smaller because the walls create a more 
uniform interference field at the model. 

EFFECTS OF COMPHESSIBnrrY ON THE CORRECTIONS 

It is shown in reference 9 that the lift due to compressible 
flow can be corrected in exactly the same manner as though 
the flow were incompressible. Thus, neither the k nor the S 
of equations (43) and (44) is a function of compressibflity. 
Since the arguments used in reference 9 are based on sub- 
sonic linearized compressible flow, it may be expected that 
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corrections can be made for the critical or slightly higher 
blach numbers. 

APPUCABIUTT OF THE THEORY 

Several of the dififerences between the idealized and the 
actual problem are those due to viscosity, which causes a 
mixing region in the neighborhood of the slot rather than 
the assumed constant pressure surface which divides the high- 
velocity tunnel air and the stagnant tank air. Since the 
mixing region involves various complicated phenomena such 
as turbulence, velocity gradients, separation at the outer 
slot edges, and differences in boundary conditions in different 
dots depending upon the direction of flow through the slots, 
it is very likely that its effects on the lift interference will 
have to be determined from analysis of experimental data 
(for example, ref. 10). 

Another effect of viscosity which becomes important if 
the slots are narrow and deep is the friction on the air as it 
flows through the slots. Since the friction reduces the 
amoimt of flow through the slots, they will effectively be- 
come narrower, so that the quality factor may be expected 
to become larger. 

The possibility also exists that, in an actual tunnel, the 
pressures may not be equal in all the slots and the slot 
pressures may not be equal to the pressure in the tank. It 
would be expected that if the difference between the pres- 
sure in the slot and that in the tank causes more flow through 
the slot than occurs in the ideal tunnel, the slot will effec- 
tively be wider and if the difference decreases the flow, the 
slot win effectively be narrower. Thus, the actual quality 
factor will depend upon the effective slot widths as deter- 
mined by the pressme increments. 

Other differences occur because a practical timnel cannot 
be constructed like the ideal tunnel. These differences in- 
volve finite slot lengths, variable-width slots, and lips on the 
slot edges. In considering the effects of finite slot lengths, 
it seems reasonable to assume that those effects should be no 
more serious than the effects of a finite-length open tunnel. 
In reference 11, it is shown that the theoretical lift corrections 
for an infinitely long open tunnel are adequate, provided the 
model is located a distance of at least one-half the tuimel 
height from the entrance and exit regions. Therefore, if the 
slot configuration is such that the width is constant over a 
section whose length is at least equal to tunnel height, the 
theoretical corrections shoifld be adequate even though the 
slot width may vary considerably outside that region. The 
effect of hps on the slots may be shown qualitatively by 
comparing the pressure gradients of the two types of slot 
configmrations, that is, with and without lips. Since the lips 
confine the flow, the pressure gradient in that configuration 


is less steep through the slot than if there were no hps. 
Therefore, the velocity may be expected to be lower, so that 
the effective width of the slots wfll be reduced and the 
quahty factor thereby increased. 

CONCLUSIONS 

An analysis of the equations which represent the inter- 
ference on the trailing-vortex system of a imiformly loaded 
wing, due to wind-tvmnel walls with mixed open and closed 
boundaries, has shown that: 

1. Slot openings of the order of 7 percent of the tmmel 
periphery for four evenly spaced slots of equal length, and 
less for larger numbers of slots, are required to reduce the 
interference on a lifting model to zero. The zero-inter- 
ference quahty factor, which is defined as the ratio of the 
interference of a slotted tunnel to the interference of a closed 
tunnel, is critical with reg>ect to the percentage of waU open- 
ing, inasmuch as a small change in waff opening wih cause 
an appreciable change in the quahty factor when its value is 
near zero. 

2. The t unn els which contain two symmetrically located 
slots showed quite different values of the interference for 
different slot locations. The differences in interference 
became smaller as the nmnber of slots increased. 

3. In the tunnels examined, a region is noted to exist about 
the center of the tunnel in which the ratio of the slotted- 
tunnel interference to the closed-tunnel interference was 
fairly constant, so that in order to obtain the corrections for 
the effects of a slotted timnel it is necessary only to multiply 
the dosed-tunnel interference by a constant. 

4. The region in which the ratio of the slotted-tunnel inter- 
ference to the dosed-tunnel interference is reasonably uniform 
was found to be larger for the tunnels containing 8 or 12 slots 
than for those containing 2 slots. 

5. An examination of tunnels containing a single slot 
showed that this slot produced a rolling moment or a cross 
flow on the modd. Either or both of these phenomena may 
be expected for timnels containing a larger odd number of 
evenly spaced, equal-width slots. These effects should, how- 
ever, decrease as the nmnber of slots increases. Similar 
interference may also be expected for any case in which the 
slots are asymmetrically arranged with respect to the model 
axes. 

6. An analysis of the effects of compressibility shows that 
the quality factor is relatively unaffected by compressibility 
throughout the subsonic region. 

Langley AlBronautical Laboratory, 

National Advisory Committbb for Aeronautics, 
Langley Field, Va., February S, 196S. 
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